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ABSTRACT 

We quantify the effect of supernova Type la peculiar velocities on the derivation of cosmological 
parameters. The published distant and local la SNe used for the Supernova Legacy Survey first-year 
cosmology report form the sample for this study. While previous work has assumed that the local SNe 
are at rest in the CMB frame (the No Flow assumption), we test this assumption by applying peculiar 
velocity corrections to the local SNe using three different flow models. The models are based on the 
IRAS PSCz galaxy redshift survey, have varying (3 — f^ 6 /b, and reproduce the Local Group motion 
in the CMB frame. These datasets are then fit for w, fl m , and 17a using flatness or ACDM and a 
BAO prior. The \ 2 statistic is used to examine the effect of the velocity corrections on the quality 
of the fits. The most favored model is the (3 = 0.5 model, which produces a fit significantly better 
than the No Flow assumption, consistent with previous peculiar velocity studies. By comparing the 
No Flow assumption with the favored models we derive the largest potential systematic error in w 
caused by ignoring peculiar velocities to be Aw = +0.04. For SIa, the potential error is Aft a = —0.04 
and for fi m , the potential error is Af2 m < +0.01. The favored flow model (/3 = 0.5) produces the 
following cosmological parameters: w = — 1.08l ; 8) ^™ = ^-^-0.02 assuming a flat cosmology, and 
O a = 0.80±g;gf and Q m = 0.27^°:^ for a to = -1 (ACDM) cosmology 



Subject headings: cosmology: large-scale structure of the universe 
- supernovae: general 



galaxies: distances and redshifts 



1. INTRODUCTION 

Dark Energy has challenged our knowledge of funda- 
mental physics since the direct evidence for i ts existence 
was d iscovered using Type la supernovae (jRiess et al.1 
119981 iPerlmutter et all I1999D . Because there are cur- 
rently no compelling theoretical explanations for Dark 
Energy, the correct emphasis , as pointed out by th e Dark 
Energy Task Force fDETF. lAlbrecht et al.ll2006fl . is on 
refining our observations of the accelerated expansion of 
the u niverse. Recommend ation V from the DETF Re- 
port (jAlbrecht et al.ll2006l ) calls for an exploration of the 
systematic effects that could impair the needed observa- 
tional refinements. 

A couple of recent studies (|Hui fc Greend 120061 : 
ICoorav fc Cald well 2006) point out that the redshift 
lever arm needed to accurately measure the universal 
expansion requires the use of a local sample, but that 
coherent large-scale local (z < 0.2) peculiar velocities 
add additional uncertainty to the Hubble diagram and 
hence to the derived cos mological paramete rs. 

Cu r rent analyses (e.g., lAstier et"aT]|2006t iRiess et all 
120071 iWood-Vasev et all [20071) of the cosmological pa- 
rameters do not attempt to correct for th e effect of local 
peculiar velocities. As briefly noted by iHui fc Greend 



(|2006h and [C oorav fc Caldwelll (|2006h . it is possible 
to use local data to measure the local velocity field 
and hence limit the impact on the derived cosmo- 
logical parameters. Measurements of the local ve- 
locity field have improved to the point whe re there 
is consistency among surveys and methods (|Hudsonl 
20031: iHudson et all 2001 iRadburn-Smith et all 12004 
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iHaugboelle et al.l (|2006h : 
IWatkins fc Feldmanl (|2007h and others. Their results 
demonstrate that the local flows derived from SNe are 
in agreement with those derived from other distance in- 
dicators, such as the Tully-Fisher relation and the Funda- 
mental Plane. Our aim is to use the current knowledge 
of the local peculiar motions to correct local SNe and, 
together with a homogeneous set of distant SNe, fit for 
cosmological parameters and measure the effect of the 
corrections on the cosmological fits. 

To produce this measurement, we analyze the lo- 
cal and distant SN la sample used in the first-year 
cosmolo gy results fr om the Supernova Legacy Survey 
(SNLS, lAstier et all I2006L hence A06). This sam- 
ple is composed of 44 local SNe (A06, Table 8: 
Hamuv et al. 1996; Riess et al. 1999; Krisciunas et al.l 
Strolger et al.l 



Jha 2002; Stro lger et all 120021 : lAltavilla et all 
Krisciunas et al.l l2004alrbf ) and 71 distant SNe 
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(A06, Table 9). The distant SNe are the largest homo- 
geneous set currently in the literature. The local sample 
span the redshift range 0.015 < z < 0.125 and were 
selected to have good lightcurve sampling (A06, § 5.2). 
Using three different models encompassing the range of 
plausible local large-scale flow, we assign and correct for 
the peculiar velocity of each local SN. We then re-fit the 
entire sample for w, Cl m , and fl\ to assess the system- 
atics due to the peculiar velocity field, and to asses the 
change in the quality of the resulting fits. 

2. PECULIAR VELOCITY MODELS 

Peculiar velocities, v, arise due to inhomogeneities in 
the mass density and hence in the expansion. Their ef- 
fect is to perturb the observed redshifts from their cos- 
mological values: czcmb = cz + v • r, where cz is the 
cosmological redshift the SN would have in the absence 
of peculiar velocities. With the advent of all-sky galaxy 
redshift surveys, it is possible to predict peculiar veloc- 
ities from the galaxy distribution provided one knows 
= / (fi) /b, where b is a linear biasing parameter relating 
fluctuations in the galaxy density, 8, to fluctuations in the 
mass density. The peculiar velocity in the CMB frame is 
then given by linear perturbation the ory dPeebles Il980f) 
applied to th e density field (see, e.g. lYahil et~ 119911 : 
lHudson|[l99l : 
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In this Letter, we use the den sity field of IRAS PSCz 
galaxies (Branchini ct al. 1999), which extends to a 
depth R m ax = 20000 km s _1 . Contributions to the pe- 
culiar velocity arising from masses on scales larger than 
Rmax are modeled by a simple residual dipolc, V. Thus, 
given a density field, the parameters f3 and V describe 
the velocity field within R m ax- For galaxies with dis- 
tances greater than R max , the first term above is set to 
zero. 

The predicted peculiar velocities from the PSCz den- 
sity field are subject to two sources of uncertainty: the 
noisiness of the predictions due to the sparsely-sampled 
density field, and the inapplicability of linear perturba- 
tion theory on small scales. Typically these uncertainties 
are accounted for by adding an additional "thermal" dis- 
persion, which is assumed to be Gaussian. From a care- 
ful analysis of pred i cted a nd observed peculiar velocities, 
IWillick k, Straussl (|1998( ) estimated these uncertainties 
to be ~ 100 km s" 1 , albeit with a dependence on den- 
sity. iRadburn-Smith et al.l (|2004f ) found reasonable \ 2 
values if 150 km s _1 was assumed in the field, with an 
extra contribution to the small-scale dispersion added in 
quadrature for SNe in clusters. Here we adopt a thermal 
dispersion of 150 km s _1 . 

For this study, we explore the results of three different 
models of large-scale flows and compare them to a case 
where no flow model is used. These models have been 
chosen to span the range of flow models permitted by 
peculiar velocity data, and all of these models reproduce 
the observed ~ 600 km s _1 motion of the Local Group 
with respect to the CMB. The first model assumes a 
pure bulk flow (model PBF, hence (3 = 0) with V having 
vector components (57, —540, 314) km s _1 in Galactic 
Cartesian coordinates. The second model assumes (3 = 
0.5 (model B05), with a dipole vector of (70,-194, 0) km 



s _1 . The third model adopts (3 = 0.7 (model B07) which 
requires no residual dipole. We compare these models 
to the no-correction scenario adopted by A06 and others 
with P = 0, V = which we call the " No Flow" or NF 
scena rio. Note that a recent comparison (jPike k Hudsonl 
2005) of results from IRAS predictions versus peculiar 
velocity data yields a mean value fit with (3 = 0.50 ±0.02 
(stat), so the B05 model is strongly favored over the NF 
scenario by independent peculiar velocity analyses. 
3. COSMOLOGICAL FITS 

Prior to the fitting procedure, the peculiar veloci- 
ties for each model are used to correct the local SNe 
(using a variation of iHui k Greend I2006L equations 11 
and 13). We then fit our corrected SN data in two 
ways using a y 2 -gridding cosmology fitter 1 (also used by 
IWood-Vasev et al.1 120071 1 . The first fit uses a flat cos- 
mology (O = 1) with the equation of state parameter w 
and Q m as free parameters. The second fit assumes a 
AC DM (w = —1) cosmology with J7a and fi m as free 
parameters. We used the same intrinsic SN photometric 
scatter (crj n t = 0.13 mag, A06) for every fit. The result- 
ing x 2 probability surfaces for both fits are then further 
constr ained using the BAO result from lEisenstei n et al.l 
(2005). The final derived cosmological parameters are 
then used to calculate the y 2 for each fit (see A06, § 5.4). 

The fitting procedure employed here differs in imple- 
mentation from that used in A06. Three additional pa- 
rameters, often called nuisance parameters, must be fit 
along with the two cosmological parameters. These pa- 
rameters are the constant of proportionality for the SN 
lightcurve shape, a s , the correction for the SN observed 
color, p c , and a SN brightness normalization, M.. We dis- 
tinguish (3 C from the /3 used to describe the flow models 
above. A06 used analytic marginalization of the nuisance 
parameters a s and (3 C in their fits. Here these parame- 
ters are fully gridded like the cosmological parameters. 
This avoids a bias in the nuisance parameters that re- 
sults because, in the analytic method, their values must 
be held fixed to compute the errors. The result is that 
our fits using the NF scenario produces slightly different 
cosmological parameters than quoted in A06. 

4. RESULTS 

The results of the cosmological fits for each model are 
listed in Table [T] and plotted in Figure Q] and Figure [2j 
They demonstrate two effects of the peculiar velocity 
corrections: a change in the values of the cosmological 
parameters, and a change in the quality of the fits as 
measured by the \ 2 statistic. 

We expect, if a given model is correct, to improve the 
fitting since our corrected data should more closely re- 
semble the homogeneous universe described by a few cos- 
mological parameters. The % 2 of the fits for each flow 
model can be compared to the y 2 for the NF scenario 
(shown by the dashed line in the figures) as a test of this 
hypothesis. Using Ay 2 = —2\nL/L^p, where L is the 
likelihood, we find that the pure bulk flow is over 10 3 
times less likely than the NF scenario, while the B05 and 
B07 models are 13.5 and 8.6 times more likely, respec- 
tively. 

We also use these data to assess the systematic er- 
rors made in the parameters if no peculiar velocities are 

1 http: //qold. astro .utoronto . ca/conley/simple_cosf itter/ 
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TABLE 1 

Peculiar Velocity Model Parameters and Results 
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results quoted in A06 marginalizing analytically over a B and /3 C (see §0 b best fit value from lPike fc Hud son (2005) 

accounted for. The largest of these are obtained by com- 
paring the B07 model with the NF scenario. This com- 
parison yields Awbo7 = +0.040 and AfiA,B07 = —0.039. 
The same comparison for the B05 model, which is only 
slightly preferred by the x 2 statistic over model B07, pro- 
duces Aw B0 5 = +0.02 7 and AQ^ B 05 = -0.026. The 
systematic offsets for f2 m are all 0.004 or less, demon- 
strating the insensitivity of this parameter to peculiar 
velocities. This is due to the BAO prior which is insensi- 
tive to local flow and provides a much stronger constraint 
for O m than for w or J7a (see A06, Figures 5 and 6). 

5. DISCUSSION AND SUMMARY 

The systematic effect of different flow models is at the 
level of ±0.04 in w. This is smaller than the present level 
of random error in w, which is largely due to the small 
numbers of high- and low-redshift SNe. However, com- 
pared to other systematics discussed in A06, which total 
Aw = ±0.054, the systematic effect of large-scale flows 
is important. IWood-Vasev et alj (|2007l Table 5) list 16 
sources of systematic error which total Aw = ±0.13. 
Aside from three method-dependent systematics and the 
photometric zero-point error, they are all smaller than 
the flow systematic. As the number of SNe continues to 
increase, and understanding of other systematics (e.g. 
photometric zero-points) improves, it is possible that 
large-scale flows will become one of the dominant sources 
of systematic uncertainty. 

The peculiar velocities of SN host galaxies arise from 
large-scale structures over a range of scales. The compo- 
nent arising from small-scale, local structure is the least 
important: it is essentially a random variable which is 
reduced by y/~N. More problematic is the large-scale co- 
herent component. Such a large-scale component can 
take several forms: an overdensity or underdensity; a 
large-scale dipole, or "bulk" flow. 

The existence of a large-scale, but local (< 7400 
km s^ 1 ) un derdensity, or "Hub ble Bubble" was first 
discus sed by IZehavi et al.1 (|1998l ). Recently Uha et ail 
( 2006) have re-enforced this claim with a larger SN data 
set: they find that the difference in the Hubble constant 
inside the Bubble and outside is AH/H = 6.5 ± 1.8%. 
If correct, this could have a d ramatic effect on the de- 
rived cosmological parameters (|Jha et al.l 120061 Fig 17), 
especially for those studies that extend their local sample 
down below z < 0.01 5. However, the "Hubb le Bubble" 
was not confirmed bv lGiovanelli et ail (|1999l ) who found 
AH/H = 1.0 ± 2.2% using the Tully-Fish er (TF) pe- 
culiar velocities, nor bv lHudson et al.l (|2004D who found 



Fig. 1. — Parameter values for the w, Q m fit (SI = 1 + BAO prior) 
for each of the four peculiar velocity models in Table [T] The values 
for the NF scenario are indicted by the dashed lines. The largest 
systematic error in ui compared with the NF fit is +0.040 for the 
B07 model, which demonstrates the amplitude of the systematic 
error if peculiar velocity is not accounted for. The offsets for Q m 
are all within ±0.003 snowing that this parameter is not sensitive 
to the peculiar velocity corrections due to the BAO prior. The x 2 
of the fits improve when using the two j3 models (B05, B07), while 
the PBF model provides a significantly worse fit. 
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Fig. 2. — Parameter values for the Cai Sim At (w = — 1 + BAO 
prior) for each of the four peculiar velocity models as in Figure [T] 
Again, comparing the NF fits to the B07 model produces the largest 
systematic in $1a of —0.039. We also find fl m insensitive to the 
corrections, having all offsets within ±0.004. The \ 2 values show 
the same pattern as in Figure [l] favoring the j3 models over no 
correction (NF), and over pure bulk flow. 
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AH/H = 2.3 ± 1.9% using the Fundamental Plane (FP) 
distances. 

According to equation [TJ a mean underdensity of IRAS 
galaxies of order ~ 40% within 7400 km s^ 1 would 
be needed to gen erate the "Hubble Bubble" quoted by 
Uha et all (|2006t ). However, we find that the IRAS 
PSCz density field of iBranchini et all ((1999) is not un- 
derdense in this distance range; instead it is mildly over- 
dense (by a few perce nt) within 7400 km s" 1 (see also 
IBranchini et al.l I1999L Figure 2). As a further cross- 
check, when we refit the lJha et al.l ()2006f ) data after hav- 
ing subtracted the predict ions of the B05 fl ow model, the 
"Bubble" remains in the Uha et all (|2006t ) data. Thus, 
the Jha et al "Bubble" cannot be explained by local 
structure, unless that structure is not traced by IRAS 
galaxies. Moreover, when we analyze t he 99 SNe within 
15000 km s" 1 from lTonrv et alj (|2003[ ) in the same way, 
we find no evidence of a significant "Hubble Bubble" 
(AH/H = 1.5 ± 2.0%), in ag reement with the r esults 
from TF and FP surveys The lTonry et all 1)20031 ) sam- 
ple and that of lJha et all (|2006[ ) have 67 SNe m common. 
The high degree of overlap suggests that the difference 
lies in the different methods for converting the photom- 
etry into SN distance moduli. 

2 http : //sdssdp47 . f nal . gov/sdsssn/sdsssn . html 



A local large-scale flow can also introduce systematic 
errors if the low-z sample is biased in its sky coverage: 
in this case, an uncorrected dipole term can corrupt the 
monopole term, which then biases the cosmological pa- 
rameters. For the large-scale flow directions considered 
here, this does not appear to affect the A06 sample: we 
note that the PBF-corrected case has similar cosmolog- 
ical parameters to the "No Flow" case. However, if co- 
herent flows exist on large scales, this may affect surveys 
with unbalanced sky coverage, such as the SN Factory 
(lAldering et al.ll2002l) or the SDSS SN survey 2 . 

The most promising approach to treating the effect 
of large-scale flows is a more sophisticated version of 
the analysis presented here: combine low-redshift SNe 
with other low-redshift pec uliar velocity tracers , such as 
Tully-Fisher SFI++ survey dMasters et al.ll2006f) and th e 
NOAO Fundamental Plane Survey ([Smith et al.ll2004h . 
and use these data to constrain the parameters of the 
flow model (f3 and the residual large-scale flow V) di- 
rectly. One can then marginalize over the parameters of 
the flow model while fitting the cosmological parameters 
to the low- and high-z SNe. 
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